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PREFACE 


A brief sumery of the research carried out under the NASA Grant 
NSG 3212 during the six nonth period ending June 30, 1981, is described 
in this report. The authorship is quoted in alphabetical order. 

T. R. Govindan is responsible for end wall flow analysis, A. Pandya has 
carried out the leakage flow neasureaent, and M. Pouagare carried out the 
hub wall flow measureaent . 

An M.S. thesis by A. Pandya is in preparation and will be sent out 
before the next reporting period. 
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NOMENCIATURE 

c 

chord langth 

^static 

(f-P . )/D 

reference 

^total 

(Pc -Pt * )/D 

5 '^reference 

^Ppltch 

(P3-P4>/D 

Cp 

'^yaw 

(P^-P2>/D 

d 

diameter of the probe head 

D 

P5-P 

e 

Internal energy 

h 

height of the probe from the wall 

P 

static pressure 

P 

(P^ + ?2 **3 ^4^^^ 


i • 1....5 pressures indicated by five-hole probe tubes (Fig. 11 ) 


inlet total relative velocity 

r,0,z 

cylindrical coordinate system 

R 

r/rt 

R 

e 

Reynolds number 

w,v,u 

velocities in radial* tangential» and axial directions* respectively 

a 

absolute flow angle measured from the axial direction 

y 

ratio of specific heats 

u 

viscosity 

c.n.r 

transformed coordinates in the streamvise* normal* and radial 
directions* respectively 

1 p 

density 

Q 

% 

angular velocity 

k 

t _ 




Subtcrlpf 

t,h values at the tip and hub, reapactivaly 

m maxlwim valua 

r,6, z values in tha radial, tangantial, and axial diractlocj, 

raspactivaly 

T 


total 


1. SOLUTION OF THE FLOW IN THE ROTOR END WALL REGION, 

INCLUDING THE EFFECTS OF TIP-CLEARANCE FLOW 

1.1 Introduction 

The end-wall boundary layer in a rotor blade row le one of the more 
conpllcated aspects of the flow in turbomachinery. The flow is three- 
dimensional and both the radial and tangential gradients of flow variables 
are liig>ortant. However, since this is a flow with a dominant flow direction, 
the streamwlse diffusion terms in the Navler-Stokes equations can be 
neglected in comparison to the diffusion in the radial and tangential 
directions. This renders the equations parabolic and amenable to solution 
by a marching procedure. 

An important part of the endvall flow phenomenon is the effect of blade 
tip clearance. Any prediction scheme of the flow in this region must take 
this effect into consideration. Figure 1 shows the flow domain of interest 
and its relation to flow in the rotor. It is immediately noticed that the 
effect of tip clearance on the solution of the flow in the domain is through 
the boundary conditions prescribed along AB and FE (Figure 1). The choice 
of AB and FE as boundaries of the domain is dictated by the nature of the 
flow and convenience, since it is obvious that these lines are neither 
physical boundaries nor are they lines of periodicity. However, it would 
seem that the flow in the region A'B'BA can be described by simple flow 
equations based on the pressure drop across the pressure and suction surfaces 
of the blade and simple modelling of the viscous terms, for the flow normal 
to the streamwise direction. Limited experimental evidence suggests that 
this is the dominant flow direction in the region. In such a situation. 
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the peraaettr that couplaa the aolutiona to the end^wall and tip^clearance 
flowB la the preaeure drop. Thia ia very convenient from the point of view 
of the aolution acheae (to be deacrlbed in later aectlona) for the end wall 
flow where a global preaeure Iteration ia uaed to update the preaaure field. 
For every one of theae iteratlona, the tip clearance flow can be updated 
using the newly calculated preaaure diatrlbutlou, and the boundary conditiona 
along AB and FE modified before proceeding to the next preaaure iteration. 

Conaequently, the firat and important part of the aolution procedure 
would be the aolution of the flow in the region ACDF. The eqtiationa and 
the method of solution are described in the next: section. 


1.2 The Equations 

The Navier-Stokes equations in a rotating cylindrical coordinate system 
can be written in the form 


3E . 3F . 3G 4. 4. il 4. Cl 


( 1 ) 


The terms Chat constitvte the vectors I,F,G,C,P,Q,T,S are described in 
Appendix I. The choice of a cylindrical coordinate instead of a Cartesian 
system was dictated by a simpler grid generation scheme in the former case. 
Two dimensional grids can be stacked in the radial direction in preference 
to the complication of generating three dimensional grids. The equations 
are transformed to a body fitted coordinate system and the flow domain 
mapped to a sector of a cylinder in which the computations are carried out. 
This transformation can be described as 


C S 5(z,9) 
n = n(z,0) 


■■ f 


r = r 
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where C the coordinece on the body surface in the streanwise direction 
and n nonsal to it on the cylindrical surface of radius r. 

The diffusion terns in the ^ direction can be dropped to nake the 
equations parabolic in this direction and anenable to solutions by a marching 
procedure. It can easily be shown that the form of the transformed equa- 
tion can be kept similar to equation (1) 


where 


and 


AAA 



A A 


AriP +12 

Re‘ 3r 3n 


A 

+ s] 


E • CjE + CqF 
? - 5,E + ZqT 

A 

G i G 

A 

c = c 

A 

P i P 

A 

s = s 


( 2 ) 


remain unchanged in the transformation. The terms that constitute the above 
vectors are described in Appendix 1. Equation (2) describes the flow in 
the computational domain and is the basis for the numerical solution procedure 
described next . 


1 . 3 Numerical Scheme 

The numerical scheme to be described is based on the Linearized Block 
Implicit Method of Briley and McDonald [1]. This procedure was chosen in 
favor of the delta scheme of Beam and Warming in that it provided more 


4 


versatility in configuring the schene for the fora of the equation described 
by equation (2). Further, the scheme is consistent in its intermediate 
steps allowing physical boundary conditions to be applied. Unlike in a 
time-marching method where only the final steady solution is important, 
this scheme must march along a spatial diraction and every step describes 
the development of the flow in the streamwise direction. 

For the purpose of describing the numerical algorithm, equation (2) 
can be written in a simpler form as 


where 


- D(q) + S(q) 


H(q) 

D(q) 

S(q) 


E(q) 

-Lrii + i2i _ ii . 

Re 3r 3n 3r| 3r 

S(q) - C(q) 


(3) 


D(q) contains the spatial derivative terms and is mult: -dlmensionax, q is 

A 

vector of dependent variables, and S(q) is a column vector dependent on q. 

C is the time-like coordinate and (n,r) are the spatial coordinates. We can 
write a second order accurate finite difference equivalent to (2) in the 
form 


- h")/AC - 6[D(q) + S(q)]"'^^ + (1 - B)[D(q) + S(q)]" (4) 


when B • 1/2 we have the Crank-Nicolson scheme. The system of equations 
described by equation (4) is still non-linear and needs to be linearized as 


follows : 
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^ ♦ o(«)2 

- s” + - q") + o«t)^ 

B-+1 . o” + - q") ♦ 0(M)* 

The llneerized schene it then given by, 

, n+1 n-,.. «n . S*» or®D j. 35% . n+1 n. 

{-§^)(q - q )/A5 - 0 + s + e{-^ ^ 

or 

(A + • AUD(q)" + S“l 

where 



and 


L 



,n+l 


rt+1 n 

q - q 


( 5 ) 


Expanding the operators into their nulti-diaentional fora 
lA + + D^)q^ + s“l 

Factorising the scheme based on the Douglas-Gunn split, we have 
lA + + D^)q" + s"] 

(A + At • A4>* 


( 6 ) 
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2 

The difference schene is of order (AO and is consistent at its intemediate 
step. It should be noted that in writing the above factorized scheme, mixed 
derivative terms occurring in the Jacobian matrices of L need to be 
explicitly differenced to keep the system of algebraic equations describing 
the discretization of equations (6) block tridiagonal. As has been shown by 
Bean and Warning [2], this does not impaiw the stability or accuracy of the 
scheme. 

It must be noted that the general solution to the equation ( 3) would 
allow an exponential term that would grow if an eigenvalue of the Jacobian 
matrix of H, 3H/3q, is negative. As has been shown by Schiff and Steger [3], 
this is indeed true for subsonic flow. Details of the analysis are in [3] 
and are not included here for brevity. This condition is to be expected 
from Che known elliptic nature of the pressure field. Schiff and Steger 
show that the eigenvalues can all be made positive by removing the pressure 
dependency frc.v the Jacobian matrix 3H/3q. The eigenvalues are positive 
provided Che screamwlse velocity is positive, restricting the method to 
unseparated flows. This, again, is an expected result. 

To remove the dependency of the Jacobian matrix of H on the pressure 
requires the specification that p ^ p(q) i.e., Che pressure field must be 
specified in the marching equations. A global pressure iteration must be 
implemented to converge the solution. To begin with, the boundary layer 
approximation, 3p/3r ■ 0 is used and the pressure field equated to that in 
the outer inviscid flow. On marching the solution a new pressure field is 
obtained that is used to update the specified pressure field. This is the 
technique of global pressure iteiat* t . 

As mentioned in a previous section, the global pressure iteration can 
also be used to update the tip-clearance flow solution. The new pressure 
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field is used to predict a new tip clearance flow solution that nodifies the 
boundary conditions to the end-wall boundary layer equations. 

1 . 4 Status of the Program 

A conputcr program has been written to implement the solution procedure 
Just described. Initial debugging of the code is complete. The present 
version is configured for a maximum 50 x 50 grid size, runs on The Pennsylvania 
State University's IBM 3033/370 processor using the Fortran H-extended 
compiler. Run times per marching step are of the order of two seconds for 
a 20 X 20 grid. This run time includes I/O time used for disk I/O operations. 

Further reductions in this time should be expected in final versions 
of the code. 

The code is currently being used to predict the flow in a 90* bend 
using the experimental data of Taylor et al. [4], Preliminary results are 
encouraging. 



2. LEAKAGE FLOW MEASUREMENT AT THE 
TIP OF A COMPRESSOR ROTOR BLADE 


Detailed hot wire iiieasurements inside the clearance region of the 
Axial Flow Compressor Rotor in the Department of Aerospace Engineering was 
undertaken recently. The rotor was operated at the design flow coefficient 
as well as one off-design flow condition. The probe employed and the 
schematic of the probe set up are shown in Figures 2 and 3, respectively. 

2,1 Clearance Measurements 

Instantaneous clearance during the motion of the rotor was found using 
a fotonic sensor, and these measurements were compared with the static 
clearance for one blade. This process was repeated for four different blades 
and no appreciable difference was found between the static and the dynamic 
clearances in all the four cases. Therefore, it was assumed that the 
clearance did not vary with the rotational speed of the compressor and 
static clearance measurements were carried out for all the 21 blades at 
various axial locations. Table 1 shows the clearance in thousands of an inch 
at five axial locations. 

2 • 2 Flowfield Measurements 

Flowfield measurements inside the clearance region and the rotor exit in 
the end wall region were carried out using a stationary **v” configuration hot- 
wire probe in combination with ensamble averaging technique. 

Figure 2 shows the two sensor hot wire probe used for these measurements. 
Unlike the conventional two sensor probes (x configuration), this probe has 
both the sensors in the same plane and are at right angles tc each other. 
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The probe was introduced in the flow field through the holes made at several 
axial locations on the lucite window. Care was taken to avoid a sudden step 
(and thereby generation of eddies) at the location where the probe was 
introduced. Figure 3 shows the arrangement made to avoid the step. Figure 4 
shows the schematic of the instrumentation used for the hot wire measure^nents 
and data processing. 

The sensors were connected to two DISA 55M10 constant temperature 
anemometer units. The d.c. components of the anemometer signals were 
measured with an integral digital voltmeter. The fluctuating components of 
these signals were amplified and recorded with SAVRE IV FM signal recorder/ 
reproducer at the tape speed of 15 i.p.s. 

The r.p.m. of the rotor was monitored by a photo cell circuit which 
used a 60 slot calibrated disk mounted on the rotor shaft. The output of 
the photocell circuit was displayed on a digital counter. The photocell 
circuit also provided one sharp pulse for each revolution of the rotor. 

This probe was also recorded along with the a.c. signals to identify a specific 
blade passage. 

The instantaneous signals from the anemometer and the pulse were 
displayed on Tektronix type RM361A storage oscilloscope. 

Measurements were taken at five axial locations for tip leakage flow 
and four locations for rotor exit flow. Rotor inlet flow was measured at 
one location 1/8 in. upstream of the leading edge. The axial and radial 
locations at which the measurements were taken are listed in Table 2. At 
each axial location, the probe was traversed as near the blade as possible 
(0.013 in. away from the blade tip with minimum clearance) without damaging 
the probe. The probe was traversed at an interval .005 in. The same 
procedure was repeated for the measurements at the off design conditions. 



Table 2, Axial and Radial Locations for Flow Field Measurements 
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For the rotor exit flow measurements, in the region away from the wall, 
the probe was traversed at the interval of .025 in. In the region very near 
the wall the probe was traversed at the interval of .010 in. 

Table 3 lists the various parameters at the two operating conditions 
at which the measurements were carried out. 

2. 3 Results 

The data at z/c ■ 0.75, 0 * 0.55 has been processed using the ensemble 
average technique. All other data will be incorporated and interpreted in 
the M.S. thesis by A. Pandya. 

Figure 5 shows the tangential variation of mean axial velocity for one 
blade passage. The mean axial velocity was normalized by the local maxiraxim 
axial velocity. It can be noted from this figure that curves for various 
radial location collapse into a single curve. Attempts are presently being 
made to identify the location of gap region in this curve. Figure 6 shows 
the variation of the mean tangential velocity which is normalized with 
respect to the local maximum axial velocity. A similar trend is observed 
as in the case of mean axial velocity i.e., a sharp decrease in the mean 
tangential velocity and all the curves collapsing into one curve in the 
clearance region. Figure 7 shows the variation of alpha (ot) , the angle made 
by the total velocity with the axis of the rotor. Once again, the clearance 
region is characterized by a sharp decrease in ol. There is a large varia- 
tion in X (approximately 27 degrees) inside one blade passage. Figures 5, 
h, and 7 show a peculiar behavior of the flow in the passage, i.e., two 
peaks in each of these three quantities. This sort of a behavior may be 
due to the scraping vortex as well as leakage flow. 
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Table 3. Compressor Operating Conditions 


1 

j 

i 

Design Condition 

Off-Design Condition 

R.P.M. 

1100 

935 

Inlet Static Pressure 

-5.08 ctns of H 2 O 

-4.356 cos of HjO 

Inlet Velocity (Axial) 

28.37 m/sec 

26.28 m/sec 

Flow Coefficient 

0.55 

0.60 


Figure 8 shows the variation of maxiiaum axial velocity at various radial 

locations. The actual distance from the casing is normalized with respect 
to the clearance height. It is evident from the figure that the maximum 
axial velocity is becoming fairly constant at 1/3 of the passage height 
(2/3 from the casing). This behavior is in accordance with the expected 
profile where the velocity reaches a maximum very near the blade tip before 
reaching zero. Unfortunately, it was not possible to go nearer to the 
blade tip with the present set-up. 

Figures 9 and 10 show the variation of the turbulence intensities in 
the axial and tangential directions, respectively. These are normalized 
by the local inaximum axial velocity. Both these quantities show a behavior 
which is very similar to that of the mean quantities. A sudden decrease in 
the turbulence intensities in the clearance region is observed. Reasons for 
very peculiar and erratic behavior at the radial station 0.011 in. from the 
casing are not known at the present time. 














3. MEASUREMENTS OF THE ROTOR HUB WALL BOUNDARY LAYER 


The flow in the hub region of a rotating blade row is one of the least 
understood aspects of the flow. A better understanding of the flow phenomena 
in this region is important for the design of the blade row, reduction of 
the losses, and improvement of the performance of the compressor. Such 
measurements are also needed to check flow prediction schemes. 

This report presents flow measurements taken near the hub wall inside 
the rotor passage of the single stage axial flow compressor facility at the 
Department of Aerospace Engineering. A five-hole probe rotating with the 
blade row is used to measure the three components of the mean velocity, 
static and total pressure. The probe is shown in Figures 11 and 12. Measure- 
ments were taken at four radial (r/r^ * 0.51, 0.5136, 0.517, 0.52) and 
four axial locations (z/c ■ 0.245, 0.45, 0.66, 0.866) at a flow coefficient 
equal to 0.555. The results are shown in Figures 13 to 16. It can be seen 
from the axial component of velocity profile that the thickness of the hub 
wall boundary layer is very small. The radial location which is nearest to 
the wall is in the outer region of the boundary layer, although it corresponds 
to a distance of 0.635 cm (1/4 in.) from the wall. The first measuring 
station is at a distance of 0.635 cm from the suction surface and the last 
at a distance of 0.79 cm from the pressure surface. At the first two axial 
locations no blade boundary layer can be seen, but at z/c * 0.66 and z/c ■ 
0.866 the outer region of the blade boundary layer on the suction side can 
be clearly identified. 

The measurements show almost ze**o radial velocity at the first two 
axial stations. At the last two axi ^ locations radial component of velocity 
exists. It can be seen that W is stronger in the suction side region than 
in the pressure side region. This is due to the fact that on the suction 
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side, the radial velocity induced by the secondary vorticity and the radial 
component of velocity in the blade boundary layer are in the same direction. 
The opposite happens on the pressure side. 

The whole set of measurements were taken very near the wall, so the 
"wall vicinity" effects were not small. In order to evaluate this effect, 
the probe wai< calibrated with the probe axis parallel to a knife-edge plate 
in a jet. The probe is moved normal to the plate. The results are shown 
in Figure 17. All the raw data were corrected accordingly. 
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Appendix 1* The Flux Vectors and Jacobian Matrices 
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where 

U - + ^qV 

V - T]^U + HqV 

To obtain the flux vectors in the untrans formed coordinate system, we 
put - 1, ^0 • 0, - 0, Tig « 1. In deriving the energy equation, heat 

transfer terms from the temperature gradient have been neglected to reduce 
the complexity of the equations. Inclusion of these terms should present no 
additional difficulty to the solution procedure. The equation of state used 
to close the system of equations is 

P/p - (Y - l)e^ 

'I o 

" (e/p) - 0.5(0*" + v" + w^) 


where 
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Tht J«cobi«n Matrix of th« V«ctor E 
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*13 “ 
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*21 " 


- ^0“^ 

*22 ‘ 

+ CgV 

*23 ’ 

^0“ 




^ 2 

*31 " 
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- CgVW 
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5,u + 

*51 * 

- -<«_ 
1 

(C 

J : 

*52 

^z <* 

+ p)/P 

*53 ■ 

Ce (e 

+ P)/P 

*55 '■ 

'•z“ 

" v> 


All other elements of the matrix are zero. 



The Jecobiiin Matrix of the Vector F 
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2 2 

-n^u - ri0uv + 

2n,u + HaV - (Y - l)n,u 

HgU - (Y - Dn^v 

-(Y “ Dn^w 

n,(Y - 1) 
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2 2 

- HgV ngfji 

n^v - (Y - DHqV 

n^u + 2TigV - (Y - 

-(Y - DHqW 

-n^uw - r)gWV 

T)^W 

HgW 

n^U + HgV 

2 

(n^u + ngv)(2ii - Ye/p) 
2 

n^(Ye/P - 0 ) - (Y - 1) 
r.g(Ye/c - $“) - (Y - 1 ) 

Y(n u + riov) 

z ^ 


All other elen»ents of the matrix are zero. 


The Jacobian Matrix of the Vector G 


where 
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*34 " 
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2 
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853 * -(Y - l)vw 

854 = (Ye/p - - (Y - l)w^ 

S35 “ Yw 

’ 2 2 2 
- 0.5(Y - 1)(u + V + w ) 


All other olenents of the matrix are zero. 


The Jacobian Matrix of the Vector C 
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All other elements of the matrix are zero. 
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The Jacobian Matrix of the Vector P 


>■21 ■ \ “n'?” 
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w 5w 
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r w4w ^u5 1 Sv 

°r D ~ ^ ^ P '2 0 ” 3 ^9 p 3n 


3 3 3 3 

D , D denote operators of the fona-;r"{( , respectivelv . 

r n dr 3n 3^ ’ 
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The Jacobian Matrix of the Vector 
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Th« Jacobi^a Matrix of th« Victor ^ 
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Fig. 1 Endvall Flow 
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Figure 2. Two-Sensor Hot-Wire Probe 
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Figure 3. Probe Set-Up for Leakage Flow Measurement 
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VARIATION OF AXIAL VELOCITY 



0.00 0.80 1.00 
TANQOCrtAL 0I8TANCK 


Radial Station 

A • 0.04” from the wall 
B ■ 0.035" from the wall 
C ■ 0.040” from the wall 
D ■ 0.025” from the wall 
E ■ 0.015" from the wall 


Figure 5. Tangential Variation of 
Velocity Inalde the Gap 
Radial Locations (z/c “ 
Clearance Height 0.059" 


Axial 

at Various 
0.75, Tip 
to 0.069") 
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Figure 7. Tengcntlel Variation of Abaolute 
Angle Inalda the Gap at Varloua 
Radial Locationa (z/c ■ 0.75, 

Tip Clearance Hai^t 0.059" to 0.069") 


Radial Station ' 

A “ 0.04" fron the wall 
B - 0.035" froa the wall 
C " 0.040" froa the wall 
D - 0.025" froa the wall 
E " 0.015" froa the wall 
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TURBULENCE INTENSITY IN TANGENTIAL DIRECTION 





Figure 12. The Geometry of the Five-Hole Probe 
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BLADE SURFACES 



Figure 13. Blade to Blade Distribution of Axial, Tangential, 
.’’nd Radial Relative Velocity in the Hub Region at 
z/c - 0,247 
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Figure 14. Blade to Blade Distribution of Axial, Tangential, 
and Radial Relative Velocity in the Hub Region at 
Z'Z • 0.45 


1 



42 






• r/r^ ■ 0.51 

* r/r^ • 0.5136 
+ r/r^ • 0.517 
o r/r^ - 0.52 


1* 


t ^ 


+ t + 


‘Z» ^ 


A Gfc ^ ,> ,i> r* » :•* ^ 






^ ^ ^ ^ 

<.':UGLfJ.\L PACK LS 

tjp p(K»K <;UAi rn 


' Fij^ure 15. Blade to Blade Distribution of Axial, Tangential, 
SS ■’nd Radial Relative Velocitv in the Hub Region at 

r/c - 0.66 
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$S Figure 16, Blade to Blade distribution of Axial, Tangential, 
and Radial Relative Velocity in the Hub Region at 
2/c « 0.866 
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Flf»ure 17. "IJall Vicinity" Effect 



